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Abstract 

We discuss the production of multi-photons squeezed states induced by the time variation 
of the (Abelian) gauge coupling constant in a string cosmological context. Within a fully 
quantum mechanical approach we solve the time evolution of the mean number of produced 
photons in terms of the squeezing parameters and in terms of the gauge coupling. We com- 
pute the first (amplitude interference) and second order (intensity interference) correlation 
functions of the magnetic part of the photon background. The photons produced thanks to 
the variation of the dilaton coupling are strongly bunched for the realistic case where the 
growth of the dilaton coupling is required to explain the presence of large scale magnetic 
fields and, possibly of a Faraday rotation of the Cosmic Microwave Background. 
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The squeezed states formalism has been successfully applied to the analysis of tensor, 
scalar Jl| and rotational @] fluctuations of the metric by Grishchuk and collaborators. In 
this paper we want to address the possible application of the squeezed states formalism to 
the case of relic photons which, to the best of our knowledge, has not received particular 
attention. In the case of relic gravitons and relic phonons the analogy with quantum optics 
is certainly very inspiring. In the case of relic photons the analogy is even closer since the 
time variation of the dilaton coupling plays directly the role of the laser "pump" which is 
employed in order to produce experimentally observable squeezed states || . 

In a general relativistic context the evolution equations of the (Abelian) gauge field 
strength is invariant under conformal (Weyl) rescaling of the metric tensor. If we are in a 
conformally flat background of Friedmann-Robertson- Walker (FRW) type (written in con- 
formal time 77) 

ds 2 = G^dx^dx" = a 2 (r])[dr] 2 — dx 2 }, (1) 

the evolution of the Abelian gauge field strength implies that, in vacuum, the determinant 
of the FRW metric can be always reshuffled by appropriately rescaling the gauge fields so 
that the (rescaled) magnetic and electric field amplitudes will always obey, in the curved 
background of Eq. (HD the same Maxwell's equations they would obey in flat space-time. 
Consequently, in general relativity Abelian gauge fields cannot be amplified from their vac- 
uum fluctuations |4], [5]. On the contrary, tensor fluctuations of the metric are amplified in 
an isotropic FRW space since their evolution equation is not invariant under Weyl rescaling 
of the background metric 

A peculiar feature of cosmological models inspired by the low energy string effective action 
is that the gauge coupling is not really a constant but it evolves in time ||. Therefore, 
in string cosmology, photons |7] (as well as gravitons ||) can be produced thanks to the 
time evolution of the gauge coupling from the electromagnetic vacuum fluctuations. A 
similar amplification effect exists for the case of non- Abelian gauge fields which are, however, 
screened at high temperature. The large (but finite) value of the conductivity in the early 
Universe implies that the Abelian magnetic component survives only for momenta which 
are much shorter than the magnetic diffusivity scale at each epoch |TI], |TT|. Plasma effects 
connected with the evolution of large scale magnetic fields have been explored with particular 
attention to the electroweak scale The effective action of a generic Abelian gauge field 
in four space-time dimensions reads 

S = -~ / d*x^G- 2 F aP F a/3 , (2) 

where F a p = V[ a Ap] is the Maxwell field strength and V Q is the covariant derivative with 
respect to the string frame metric G^ IU . In Eq. (|2|) g = exp (0/2) is the (four dimensional) 
dilaton coupling. 

In the absence of a classical gauge field background, the perturbed effective Lagrangian 
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density 



z a L y 



, 1,(77) = y d 3 xC(x, 77), 



(3) 



describes the evolution of the two ( a — <g>, ©) transverse degrees of freedom defined by 
the Coulomb gauge condition A = and V • A = (the prime denotes differentiation 
with respect to conformal time). The fields A a = gA a have kinetic terms with canonical 
normalization and the time evolution 



^-V 2 A a -g(g- 1 )"A a = 



(4) 



can be derived from the Euler-Lagrange equations By functionally deriving the the action 
we get the canonically conjugated momenta ir a = A' a + (g'/g)A a leading to the Hamiltonian 
density and to the associated Hamiltonian 



7C 2 a + (diAaf ~ 1-AoJXo, 
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H(rj) = J d 3 xH(x). 



(5) 



The operators corresponding to the classical polarizations appearing in the Hamiltonian 
density 



A a (x,r)) = J d 3 k ^ 3/2 A a (k,ri) 



Ak-x 



Aa(k,r]) = -^[akM + a-kM]* 



7t a (x,r}) = J 



'k 



(2tt) 3 / 2 



Tt a (k,7])e lk - X , 



(6) 



obey canonical commutation relations and the associated creation and annihilation operators 
satisfy aj,^] = 8 a p8^\k — p). The Hamiltonian can then be written as: 



H(v) = d 3 kJ2 k (4,<A,a + al kia a- k , a + 1) 



+e(y)a_ fc , Q a fciQ + e* {g)a{ a a ] _ kf 



, e(y) 



■ 9' 



(7) 



The (two-modes) Hamiltonian contains a free part and the effect of the variation of the 
coupling constant is encoded in the (Hermitian) interaction term which is quadratic in the 
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creation and annihilation operators whose evolution equations, read, in the Heisenberg pic- 
ture 

— — = «fca^ a_ fe)Q . 8 

The general solution of the previous system of equations can be written in terms of a 
Bogoliubov-Valatin transformation 

dk, a (v) = A i *,a(»7)6fc ) a + "kAvW-^a 

4,a0?) = /4,afa)&I,a + "jU^"*." ( 9 ) 

where a& !a (0) = and a_& jQ ,(0) = &_& ia . Unitarity requires that the two complex functions 
//fc(?7) and v k {j]) are subjected to the condition \nkijf)\ 2 — \vk(j])\ 2 = 1 which also implies that 
fik{v) an d ^k(v) can De parameterized in terms of one real amplitude and two real phases 

fi k = e i9k coshr k , v k = e i(2 ^ _efc) sinhr fc , (10) 

(r is sometimes called squeezing parameter and <p k is the squeezing phase; from now on we 
will drop the subscript labeling each polarization if not strictly necessary). The total number 
of produced photons 

(Q-k%\a\{ri)a k {rj) + a j _ k a_ k \0 k 0_ k ) = 2 n k . (11) 

is expressed in terms of n k = sinh 2 rfc, i.e. the mean number of produced photon pairs in 
the mode k. Inserting Eqs. (j^),(|10D and (|TT|) into Eqs. (|j) we can derive a closed system 
involving only the n k and the related phases: 

^ = -2/(rI fc )-cos20 fc , (12) 
dr] g 

*F = "* + i7S) sm2fe (13) 

* -i+^Bin%, (14) 



di] g drik 



where f(n k ) = \Jn k {n k + 1). 

In quantum optics |T2| the coherence properties of light fields have been a subject of 
intensive investigations for nearly half a century. Magnetic fields over galactic scales have 
typical frequency of the order 10~ 14 -10™ 15 Hz which clearly fall well outside the optical range. 
Thus, the analogy with quantum optics is only technical. The same quantum optical analogy 
has been successfully exploited in particle [|13| and heavy-ions physics fl4j of pion correlations 
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(in order to measure the size of the strongly interacting region) and in the phenomenological 
analysis of hadronic multiplicity distributions. 

The interference between the amplitudes of the magnetic fields (Young interferometry 



T5|, in a quantum optical language) estimates the first order coherence of the magnetic 



background at different spatial locations making use of the two-point correlation function 
whose trace over the physical polarizations and for coincidental spatial locations is related 
to the magnetic energy density. Recall that in our gauge 

^9 ST^ „a f r. - rs Jk-x -t „-i£-#l 



Bk ^ V) = (2n)W(r,) ^ 6t J ^ k ^K« eik ' X ~ 4,^1. (15) 



By summing up over the polarizations according to 

^ = E c ?(*) e iW=(*v-^). (16) 

we get that 

Gi^v) = (0_ fe fe |A(^ry)4(f +r,77)|O fe O_ fc > (17) 



ik-r 



can be expressed, using Eqs. (g) and (|T0|) 

Qijir) = J d 3 kG tJ (k)e 
where 

+ sinh 2rfc cos 20fe] (19) 

(the vacuum contribution, occurring for r k = 0, has been consistently subtracted). The 
intercept for r = of the two-point function traced with respect to the two polarizations is 
related to the magnetic energy density 

^Pb_ _ g\vP 4 p sinh2 ^ + ginh cog (2Q) 



din a; 2n 2 

(where u = k/a is the physical frequency). The two-point function and its trace only depend 
upon rlk and upon </>/.. Since Eqs. dl2]) and ([14]) do not contain any dependence upon 9 k 
we can attempt to solve the time evolution by solving them simultaneously. In terms of the 
new variable x = kt] Eqs. (0) and ([14]) can be written as 

d(j) k d\ngdf(n k ) . 

— — = -l + — __ sin 20 fc , (21) 

ax ax anfc 

,/ " A -2/(n fc )cos20 fc , (22) 



din g 
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If \(d]ng/dx)(df(nk)/dn] C ) sin 20^1 > 1, then Eqs. (pi]) and (22) can be written as 

du k _ 2 dJ^ u ^ ^L. = _ 2m) lz^ (23) 



d In g drik ding 1 + u k 

where <j)k = arctanw^. By trivial algebra we can get a differential relation between u k and n k 
which can be exactly integrated with the result that u\ — f(n k )uk + 1 = 0. By inverting this 
last relation we obtain two different solutions with equivalent physical properties, namely 



Uk{n k 



[\(\/n h (nk + 1) ± y/n k (n h + 1) - 4)]. (24) 



If we choose the minus sign in Eq. ( 24]) we obtain that 0& ~ (m + l)vr/2, m = 0,1,2... 
with corrections of order l/n k . In the opposite case <fik ~ arctan (rIfc/2) within the same 
accuracy of the previous case(i.e. 1/njfe). By using the relation between u k and n k the 
condition \(d\ng/dx)(df(n k )/dn k ) sin 20^1 > 1 is equivalent to x<l, if, as we are assuming, 
\g'/g\ vanishes as rf 2 for rj — » ±oo and it is, piece-wise, continuous. By inserting Eq. ([24]) 
into Eq. ( |21| ) a consistent solution can be obtained, in this case, if we integrate the system 
between rjf and rji defined as the conformal times where \(d\ng/dx)(df(n k )/dn k ) sin 20^1 = 1: 

Mm) A-gJ^-g^i))^ 

+ + U), m = 0,1,2... (25) 
2 n k (ri) 

If \(d\n g /dx)(df(n k ) /dn k ) sm2(p k \ < 1 (i. e. x > 1) the consistent solution of our system is 
given by 

n k {rjf) = sinh 2 ^2 J In g(x') sm2x'dx'^J 

4>k ~ ~kr] + (pk, ifk — constant. (26) 

If the coupling constant evolves continuously between — oo and +oo with a (global) maximum 
located at some time r\ r then, for x > 1, n k ~ const.. Indeed by taking trial functions with 
bell-like shape for \g' / g\ we can show that n k oscillates around zero for large <ft k . 

Up to now our considerations were general. Let us give some specific example of our 
technique. In the low energy phase (rj < rj s ) of the pre-big-bang evolution the dilaton 
coupling is determined by the variation of the low-energy effective action [Q: 

a(rj) ~ l^r^+i, g(rj) ~ \rj\'^ /2 Tj<rj s . (27) 

During the stringy phase the average time evolution of the coupling constant can be described 
by: 

a{v) - 9(v) - \v\' P , Vs <V < Vr, (28) 
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where (3 = —{4> s — 0r)/(21nz s ) where z s = r] s /r) r . For rj > r/ r , the background is dominated 
by radiation (i.e. a{rj) ~ 77) and the coupling constant freezes to its constant value (i.e. 
4> = 4> r = const, for r\ > r) r ). Notice that g{j] r ) = exp (<f> r /2) = g r ~ 0.1-0.01. For kr] < 1 we 
have, from Eqs. ( p5|) that 

n k (r] r ) ~ \r]i/Vr\ 213 ~ | A:/^ | — 2/3 , k s < k < k r (29) 

where |?7i| ~ k~ l < \rj s \. Similarly, if > \r) s \, 

nM-lk/hl^lg^/g^r)]- 2 , k < k s . (30) 

Notice that we assumed (3 > which means that the coupling constant does not de- 



creases during the stringy phase. Due to magnetic flux conservation [ITUl O the fraction of 



electromagnetic energy stored in the mode u does not change and it is defined as 

\(uj) = — j Pb = ^— n k {r) r ) sin 2 kr), (31) 



Pjdlnu An \u r 

4 / r, \ 2 



PM = M*H?(^) (32) 



where cu r ~ a r /r] r = yj g r / '47T10 11 Hz is the maximal amplified frequency red-shifted today 
and where we assumed nk(r] r ) > 1. Notice that in the unifying notation of eq. ( |3lD the 
oscillating part occurs, for each mode, when krj > 1 but not in the opposite limit where 
0fc ~ 7r/2. The critical density constraint, implies, during the stringy phase that (3 < 2. 
If (3<2 (for instance f3 ~ 1.9) we can have that A^dec) = gf (^dec/^r) 4 ~ 2/3 ~ 10 -8 for 
^dec ~ 10 -16 Hz (for c^dec > w,). Recall that, in order to rotate the polarization plane of 
the Cosmic Microwave Background Radiation, we need, at decoupling, £?> 10~ 3 Gauss, or 
in our language, A(c<Jdec) ?^10~ 8 Similarly, at the scale of 1 Mpc (i.e. ujq ~ 10 -14 Hz) 
we can have \{ujg) > 10 -10 0. Both at the galactic and decoupling frequencies n w 3> 1 in 
the framework of this specific model and the quantum mechanical state is strongly squeezed 
(N >!)• 

The quantum degree of second order coherence is a measure of the correlation of the 
magnetic field intensities at two space-time points. The intensity fluctuations of a given 
light field are described by the Glauber correlation function [12| which is nothing but the 



quantum mechanical generalization of the correlation between the classical intensities of 
two light beams. In quantum optics one deals mainly with intensity correlations of electric 
fields. This is due to the fact that the photons of the visible part of the electromagnetic 
spectrum are detected via photo-electric effect, and, therefore, what is indeed detected is an 
electric current induced by a photon. In our case we are mainly interested in the intensity 
correlations of the magnetic part of the field and we can write the corresponding correlation 
r(r) of the intensity operators as 

(: (3-(x } rj) (3-(x + r,rj) (3 + (x + r,r/) (3 + (x,r]) :) 
(: p-(x,ri)$ + (S,rj) :)(: f3~-(x + r,r)) :) 
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In this (normal ordered) definition the field operators refer to a single polarization of the 
field, namely 



rj) = j^j~ 2 J d 3 k^d k ( V )e^, (34) 

with (3~ = + Y ■ r(f) describes correlation between intensities in the case where the 
photons are detected simultaneously in time but at different spatial locations. The statistical 
properties of the given quantum mechanical state of the field are encoded in the intercept of 
the Glauber function namely T(0). 

By using Eqs. (|9|) into Eq. ([!3|) we obtain 

r ,_* _ Jd 3 k k\is k (r])\ 2 J d 3 p p\is p (ri)\ 2 



+ 



/ d 3 k k\v k (rj) | 2 / d 3 p p\fi p (r}) \ 2 
Jd 3 k k\u k {r])\ 2 J d 3 p p\u p {r])\ 2 e l ^- £ 



Jd 3 k k\v k (r})\ 2 f d 3 p p\[ip(i])\ 



2 



Jd 3 k Jd 3 ppk u* k {r])fi* k {r])u p {r])fi p (r])e i( - k 

J d*k k\u k ( V )\ 2 J d 3 p p\fi p ( V )\ 2 1 } 

which in the limit |f| — > becomes 

rrm = 2 , Jkd 3 kfp <^pu^(ri)fj,* k (ri)u p (ri)fj, p (ri) 
{) J d 3 k k\is k (rj)\ 2 J d 3 p p\iy p (v)\ 2 ' 1 J 

In order to interpret this formula we can introduce a further simplification, namely we can 
restrict our attention to a single mode of the field. Then, fx k (r]) ^K(v)^ 3 Kk — K) and 
Vkiv) ~^ VK(r))5( 3 \k — K). Thus T(0) = 3 + |z7x| -2 - As we showed, large variations in 
g(t]) give 2 ^> 1 • In the case of a coherent state the intercept of the Glauber function 



is exactly one, namely T(0) = 1 [12|. This property is in direct correspondence with the 
Poissonian character of the statistics. In the case of a thermal state (i.e. "white light"), 
T(O) — > 2 ||12|| . In the case of squeezed relic photons for large number of particles in each 
Fourier mode T(0) — > 3. Since T(0) represents the probability of two photons arriving at 
the same location, this is referred to as photon bunching. Conversely, a field with sub- 
Poissonnian statistics will have T(0) < 1, an effect known as photon anti-bunching in the 
context of the Hanbury Brown- Twiss interferometry W7\ . 
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